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AN  INTEGRAL  INEQUALITY  WITH  APPLICATIONS 
TO  ORDER  STATISTICS 


by 

Philip  J.  Boland  and  Frank  Proschan 

ABSTRACT 

We  say  the  life  distribution  function  G  majorizes  the  life  distribution 
function  F  (written  G  ?  F)  if 

CO  00 

/  G(t)dt  £  /  F(t)dt  for  all  x  £  0 

x  x 

00  00 

and  ^  G(t)dt  =  ^  F(t)dt  <  +  ■»  . 

An  integral  inequality  is  proved  giving  sufficient  conditions  on  functions 

ip  and  <f  in  order  to  ensure  that  whenever  G^^  5  for  i=l,...,  n,  then 
00  00 

^  <Kt)4(G1(t),...,Gii(t))dt  5  £  IKt)*(F1(t),...,Fn(t))dt. 

Applications  in  reliability  theory  and  order  statistics  are  given. 


1.  Introduction. 


For  given  life  distribution  functions  F  and  G,  the  respective  survival 

functions  are  F  =  1-F  and  (f  =  1-G.  We  define  the  partial  ordering  5  on  the 

m 

class  of  life  distributions  with  finite  means  by  G  >  F  (m  for  majorization)  if 


(1.1) 


and 


(1.2) 


/  G(t)dt  £  /  F(t)dt  for  all  x  >  0 

X  X 


00  <30 

UG  =  ^  G(t)  dt  =  ^  F(t)dt  —  Wp  <  ♦  00  • 


If  X  and  Y  are  nonnegative  random  variables  with  respective  distribution 

functions  F  and  G,  then  Ross  [11]  says  "Y  is  more  variable  than  X"  (written 

Y  X  or  G  ^  F)if  (1.1)  holds.  Stoyan  [14]  equivalently  defines  Y  to  be 

"larger  in  mean  residual  life"  than  X  (written  G  §c  F  or  in  previous  publi- 

(2) 

cations  G  2  F)  if  (1.1)  holds.  Bessler  and  Veinott  [3]  use  the  terminology 
"Y  is  stochastically  larger  in  mean  than  X."  The  notation  of  Stoyan  (c  for 
convex)  is  suggested  by  the  following  characterization: 

G  Sc  F 

<=> 

00  00 
Y(t)dG(t)  ^  J  Y(t)dF(t) 

holds  for  all  increasing  (that  is  nondecreasing) 
convex  functions  Y,  provided  the  integrals  exist. 

For  life  distribution  functions  F  and  G,  G  5  F  if  and  only  if  G  |c  F 
(or  G  F)  and  G  and  F  have  equal  finite  means  (Up  = y^) .  For  distribution 

functions  with  finite  means,  the  following  useful  characterization  of  G  5  f 
(see  for  example  Ross  [11]  or  Stoyan  [14])  is  an  immediate  corollary  of 
Theorem  2.1: 

G  9  F 
00 

Y(t)dG(t)  2  ^  Y(t)dF(t) 

holds  for  all  convex  functions  Y ,  provided  the  integrals  exist. 


-  2  - 


We  note  in  particular  that  if  G  ?  F,  then 

eo  co 

°G  =  (  Ct-uG)2dG(t)  £  J  (t-uF)2dF(t)  =  o2p 


Hence  G  >  F  implies  that  the  life  distribution  represented  by  G  is  'more 
dispersed*  than  that  represented  by  F  around  their  common  mean. 

For  life  distribution  functions  F  and  G  with  a  common  mean,  G  9  F  is 
a  more  general  relationship  than  G  5  F  (F  is  star  shaped  with  respect  to  G) . 

When  F  and  G  are  continuous  life  distributions  (where  F(0)  =  G(0)  =  0,  F  and 
G  have  interval  support  and  G  is  strictly  increasing  on  its  support) ,  then 
G  }  F  if  G-1F(x)  is  star-shaped  (that  is  G  i-s  increasing  for  x  >  0) . 

If  G  J  F  and  F  and  G  have  a  common  mean,  then  ¥(x)  crosses  G(x)  once  and  from 
above  as  x:0  -*■  <»,  so  that  in  particular  G  9  F  (see  Barlow  and  Proschan  [2]). 

For  a  continuous  life  distribution  function  F  with  mean  u,  let  us  define 
G(x)  =  1-e  to  be  the  exponential  distribution  with  the  same  mean.  Then 
F  is  I FRA  (increasing  failure  rate  average)  <?=»  G  >  F,  and  F  is  HNBUE  (harmonic 
new  better  than  used  in  expectation)  <=>  G  ?  F.  See  Klefsjo  [6]  for  further 
properties  of  HNBUE  distributions. 

If  F  and  G  are  two  life  distribution  functions  with  common  mean  and 
F(x)  crosses  G(x)  once  and  from  above  as  x:0  -►  ®,  then  G  9  f,  however  the 
converse  is  clearly  not  true.  For  example  let  F  and  G  be  defined  as  follows: 


F(x) 


x  <  2 
2  S  x  <  4 
4  S  x 


G(x) 


V4 

% 

.  1 


1 


T 


1 

<?- 


0 1— — 

1 

Then  G  ?  F  and  G  'crosses' 


iF 

I 

I 


I  I 

I  I 


I  I 

t  I 

f  I 

. t  “ 

I 

l 

— t - 1 - i — 

2  3  4 

F  three  times. 


I 

4- 

5 


x  <  1 
1  s  x  <  3 
3  S  x  <  5 
5  S  x 
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A  vector  b_  =  (b^,...,bn)  majorizes  the  vector  a  =  (aj,...,an)  if 

n  n 

I  ^[il  a  I  a[il  ^or  k=2,...,n 

i=k  L1J  i=k  L1J 


and 


l  ^ 


i  aCi]  ’ 


i=l 


where  the  bj-^'s  and  a[i]'s  are  t*ie  components  of  b  and  a  respectively  in 

m 

ascending  order.  When  t>  majorizes  a  we  write  b  >  a. 

Suppose  now  that  Ij  and  a  are  n  dimensional  vectors  with  nonnegative 
m 

components  such  that  b  >  a.  If  G  and  F  are  respectively  the  distribution 

functions  for  the  uniform  distributions  on  the  components  of  b  and  a,  then 
m 

G  >  F.  This  is  our  motivation  for  using  the  letter  m  for  our  partial 
ordering  on  the  family  of  life  distribution  functions  with  finite  means. 


2 .  An  Integral  Inequality. 

The  following  theorem  is  a  variant  of  an  integral  inequality  obtained 
by  Fan  and  Lorentz  [4]. 

Theorem  2.1.  Let  #  =  [0,l]n  -*•  [0,»)  be  a  continuous  increasing  function, 

and  assume  that  for  i*l,...,n,  F^  and  Gi  are  life  distribution  functions 
m 

where  G.  >  F. . 

l  l 

a)  If  i>  is  nonnegative  decreasing,  $  is  convex  in  each  variable 
separately  and  4>  satisfies  the  following  property: 

(2.1)  *(ui+  h,  Uj+  k)  -  »(ui+  h,  Uj)  -  #(u^ ,  u^+  k)  +  4>(uif  u^)  *  0 

for  all  i  *  j,  0  £  u.  £  u.+  h  £  1,  0  £  u.  £  u.+  k  £  1 

ii  J  J 

(where  we  have  used  the  notational  simplification  of  omitting  those  arguments 
of  *  which  are  the  same  in  a  given  fontula) , 
then  providing  the  integrals  exist, 

m  gg 

(2.2)  f  «Kt)*(G1(t),...,  Gn(t))dt  £  ^(t)*(F1(t),...,  Fn(t))dt. 


b)  If  \fi  is  nonnegative  increasing,  ♦  is  concave  in  each  variable 
separately  and  *  satisfies  the  following  property: 

(2.3)  *(ui+  h,  Uj+  k)  -  *(ui+  h,  u^)-*^,  U..+  k)  ♦  u..)  £  0 

for  all  i*j,  0  £  iu  s  u^  h  £  1,  0  s  u.  s  u.+  k  i  1  , 

then  providing  the  integrals  exist 

00  00 

(2.4)  /  *(t)*(G1(t),...,G  ft))dt  *  /  ♦(t)*(F.(t),....,  F„(t))dt. 

0  n  0  1  n 

Proof :  We  prove  only  a),  the  proof  of  b)  following  in  a  similar  fashion. 

(i)  Initially  we  show  that  it  suffices  to  prove  the  result  for  the  case 
when  Fj,  G.,...,  F  ,  Gn  all  have  finite  support.  In  turn  to  establish  this 
we  show  that  if  the  inequality  is  valid  whenever  F^  and  Gj  have  finite 
support,  then  it  is  true  in  general. 

Suppose  now  that  F^,  G^,...,  Fn,  Gn  are  arbitrary  life  distributions 
m 

where  G^  >  F^  for  i  =  1 . n.  Given  e  >0,  we  can  find  S  so  that 

00 

/  *(t)*(G1(t),...,<T  (t))dt  <  e  . 

S  n 

Now  define  F^ '  and  G^ *  by 

Fi'(t)  *  F^t)  t  <  s 

0  t  2  S 

V(t)  3  G  (t)  t  <  s  S  _  S  _ 

/  F  (t)dt  -  /  G. (t)dt 

—  0  1  ft  1 

G  (S)  S  £  t  £  S  ♦  - - = - - - 

Gj  (S) 


0 


otherwise 


5  - 


(if  G^(S)  =  0,  then  both  G^  and  F^  have  finite  support).  Then  G^'  >  F^',  and 

CD  OD 

/'l»(t)*(F1(t),  ?2(t),...,Fn(t))dt  >  /  ♦(t)*(?1'(t),  F2(t) . Fn(t))dt 


*  /  *(t)*(G1'(t),  G2(t),...,Gn(t))dt 


s  /  i|.(t)*(G1(t),  G2(t),...  ,Gn(t))dt  -  £  . 

Since  e  is  arbitrary,  the  conclusion  follows. 

(ii)  It  now  remains  to  show  that 

00  00 

/  *(t)*(G1(t),....(L(t))dt  <  /  *(t)*(F1Ct),...fF  (t))dt 

pin  0in 

m 

whenever  G.^  >  F^  for  all  i  =  1,...,  n,  and  where  the  support  of  and 
Gj  c  [0,S]  for  all  i  =  1,...,  n. 

Let  e  >  0  be  given.  As  *  is  continuous,  there  exists  a  6  >  0  such 

that  whenever  u,  v  e  [0,l]n  and  ||  u-v||  =  max  |u.-v.|  <  6,  then 

i=l,...,n  1  1 

|  Hu)  -  *(v)|  <  e/2S*(0). 


There  exist  only  a  finite  number  of  points  r  in  [0,S]  where  at  least 
one  of  Fj,  Gj,...,Fn,  Gn  has  a  jump  discontinuity  with  jump  >  u/2.  Hence 
we  can  find  an  integer  N  large  enough  so  that 


» 


-  6 


Hence  for  each  i  =  we  define  the  following  simple  survival 

functions: 


( 


(j+DS/N 

jS/N 


Fi(t)dt) /s/N 


and 


Gi»(t) 


when  t  e 


Qs  (j>i)s] 
LN'  N  J 


Kj+ns/N  _  . 

G. (t) dt)  / S/N 
jS/N  1 


for  some  j  =  0,...,N-1,  and  zero  otherwise. 


m 


Note  that  G."  >  F."  for  all  i  =  l,...,n. 
11 

Moreover , 


s  s 

/  ^(t)$(F1Ct),...,FnCt))dt  -  J  <Kt)*(F^'  (t),...,F»  (t))dt 


N-l  f(j  +  DS/N 

l 


j=0 


'jS/N 


*(t)[«(F1(t),...,Fn(t))dt  -  *(F»  (t),...,F«  (t))jdt 


<  *(0)2r  sup  |«|  |  ♦  (J) 


<  e  . 
Similarly, 


s  s 

J*(t)*(G1(t),...,Gft))dt  -  fU t)«(G"  (t) . G"  (t))dt 

0  0 


<  £ 


Therefore,  it  suffices  to  prove  (2.2)  for  the  case  when  all  F^,G^  are  step 
functions  which  are  constant  on  jj^,  j  ,  j  =  0,  ...»  N-l.  Furthermore, 

without  loss  of  generality  we  may  assume  that  <|i  is  constant  on  each  interval 


of  the  form 


for  j  =  0 ,  . . . , 


N-l. 


-  7  - 


m 

(iii)  Assume  now  that  for  i  =  l,...,n  and  that  all  2n  functions  have 

support  in  [0,S)  and  are  constant  on  each  interval  for  j=0,...,N-l. 

We  also  assume  t|»  is  constant  on  each  of  these  intervals  and  use  the  notational 


simplification  i|»(j)  =  iJj 


fill 

N 


for  j=0, . . .  ,N-1 . 


Each  may  be  transformed  into  by  a  finite  succession  of 
transformations  t  of  the  following  type  (see  Hardy,  Littlewood  and  Polya  [5]) 


t  changes  the  value  v^  of  G^  on  the  interval 


~i£  (21 
_N  ’  N 


1)S 


into  v . . +  h  and 
Ji 


—  PkS  (k+1) S 

the  value  v,  .  of  G.  on  ■—  into  v,  .  -  h  where  j  <  k  and 

ki  l  [_N  N  ki  J 


0  s  v,  .  -  h  s  v,  .  Sv..  <  v . .  +  h  <  1 . 
ki  ki  ji  ji 


Letting  A  denote  the  change  in  the  integral  /  ip(t) $(G. (t) . G  (t))dt 

t  Qin 

resulting  from  such  a  transformation  t,  we  complete  the  proof  by  showing 
that  A^  2:  0.  Without  loss  of  generality  i=l,  and  hence 


\  ■  jf  h,  vj2,...,».n)  -  ♦(vjl,  vj2 .  vjn)] 


-«k)C(*(»M.vk2 . vkn)  -  *(vkl-  h,»k2 . vkn)]> 

1  h'vj2 . V  -  *<vjl-''j2 . V 


-C*(Vj1+  k.vk2 ....  >vkn)  -  *(vji'vk2’"’*vkn^ 


(since  $  is  convex  in  each  variable  separately) 


8 


*  *(k)fi«*(vjl*  h.vk2+  h2* 1  * '  ,vkn+  V  “  *(vjl’vk2+  h2’*‘-Vkn+  V 
-  *(vjl+  h,vk2+  h2,...,vkn)  ♦  *(vjl.vk2+  h2,...vkn)] 


+  C*(vjl+  h,vk2+  h2*Vk3‘--*'vkn)  -  *(vjl'vk2*  Vvk3*  * '  * 'V 


-  *(vji+  h’vk2’"*’vkn)  +  *(vjl»vk2 . vkn)]} 


2  0 

(since  $  satisfies  property  (2.1)  and  t|>  is  nonnegative). 

Here  h .  =  v . .  -  v.  .  for  i  =  2 , . . . ,n . 
i  ji  ki 


Corollary  2.2  .  Let  G  and  F  be  life  distribution  functions  with  finite 
means.  Then  G  ™  F  if  and  only  if 

a)  For  all  nonnegative  increasing  continuous  convex  i>  and  nonnegative 
decreasing  ip. 


and 

b) 


/  Kt)$(G(t))dt  <  /  iKt)$(F(t))dt 
0  0 

For  all  nonnegative  increasing  continuous  concave  4>  and  nonnegative 
increasing  ip, 


/  <|»(t)*(G(t))dt  2  /  i|»(t)*(F(t))dt, 
0  0 


provided  the  integrals  exist. 


Proof.  The  only  if  part  follows  immediately  from  Theorem  2.1  .  Assume 
now  a)  and  b)  hold.  Letting  $(u)=u  and  ^(t)  =  X(-x>+09j  (that  is  the 
characteristic  function  of  the  interval  [x,+«))  it  follows  from  b)  that 

09  00 

/  G(t)dt  2  /  F(t)dt  for  all  x  2  0.  Taking  *(t)=l,  it  follows  from  a)  that 
x  x 


-  9  - 


Corollary  2.3.  If  G  and  F  are  life  distributions  with  finite  means,  then 


„  m  _ 
G  >  F 


<=> 


(2.5)  /  f(t)dG(t)  *  /  t(t)dF(t) 

0  0 

holds  for  all  convex  functions  ¥,  provided  the  integrals  exist. 

m 

Proof.  The  if  part  of  the  result  is  immediate.  Now  suppose  G  >  F.  It 

suffices  to  prove  (2.5)  for  the  case  where  4f  has  derivative  ty  and  '¥(0)=0. 

Then  «  00 

/  nt)dG(t)  =  /  *(t)G(t)dt 
0  0 


j[*(t) 


-4>(0)  ]G(t)dt  +  Ip(0)u 


-  /l>(t)-*(0)]F(t)dt  +  4<(0)y  (by  Theorem  2.1) 

0  h 


=  /  V(t)dF(t). 

0 


Remark  2.4.  Another  approach  to  (2.5)  in  the  proof  of  Corollary  2.3  is 
as  follows.  Suppose  G  >  F.  Let  Z^  and  Zp  be  the  random  variables  with 


respective  densities  — / G(t)dt  and  —  /  F(t)dt.  Then  Z  >  Z_  (Zr  is 

FG0  0  G  F  G 

stochastically  larger  than  Zp)  and  hence  (see  for  example  Ross  [11]) 

E(*(Zg))  £  E (ip(Zp) )  for  all  increasing  \f>.  But 

00  00 
/  *(t)G(t)dt  =  E(«(Zg))  *  E(ij»(Zp) )  =  /  i|;(t)F(t)dt. 


1  f  nf 


St 


3.  Applications. 

Theorem  3.1.  Let  . »Xn,Y^ , . . . ,Yn  be  independent  nonnegative  random 

variables  where  and  Y^  G^  for  i=l,...,n,  and  let  X^-j,...,  X^ 

and  ^[j]» • • • »^[n]  respectively  the  X  (Y)  observations  in  increasing  order. 
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m 

Assume  that  >  F^  for  i=l,...,n.  Then 


a)  /PCY[n]+...  +  Y[k]>t]dt  £  /  P[Xrn1+...+  Xrvi>t]dt 


[n] 


v[k] 


for  all  x  2  0  and  k  *  l,2,...,n. 
b>  <EVm . EYtn]>  5  <EXm . E1W 


Proof .  b)  follows  immediately  from  a).  In  what  follows  £  =  (ej,...,e  ) 
will  denote  any  vector  whose  components  are  zeroes  or  ones.  For  i  =  1,... 


define  4>^:[0,1]  -*■  [0,+°>)  by 


1-e 


e,  e  1-e, 

=  l  u,  ...  u,  n(l-u,)  ...(1-uJ  n 


1 

e,e.+. . .+£  2  n-i+1 

-  i  n 


n 


n' 


We  note  that  EX,-.-,  =  /  (F,  (t) . F  (t))dt  for  i  =  l,...,n.  Now  for 

n  J  q  i  i  n 

k  =  1 , . . .  ,n  we  define 


n 


*k(V--*'V  =  J  ♦i(ul»***»un) 

i=k 


n 

=  I  l 

i=k  .  .  , 

£  n-i+1 


e.  £  1-e,  1-E„ 

i  n  /  <  %  i  /  <  \  n 

ux  •••  un  .  ..(l-un) 


n  e,  en  1-e,  1-e 

=  £min(j  ,n-k+l)  £  ^  ...Uj,  (1-u^  ...(1-uJ 

j=l  e,e1+...+en»j 


Since  /  PCX[n]+...+  X[k]  >  t]dt  =  /  •fc(F1(t) , . . . ,Fn(t))dt, 

it  suffices  by  Theorem  2.1  b  to  show  that  each  #k  satisfies  (2.3)  and 

is  concave  increasing  in  each  variable  separately. 

n-k  E2  e  1-e 

...u„  (1-u,] 

'1  "  j=0  e1,e2*..;+en- j 

where  e^  represents  an  n-1  component  vector  of  zeroes  and  ones 


e2  e  1-e-  1 

u2  ...un  (l-u2)  ....(l-un) 


n,  we 


i 

A 
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As  , . . .  ,un)  is  symmetric  in  u^, _ ,un>  it  follows  that  4^  is  30 

increasing  function  linear  (and  hence  concave)  in  each  variable  separately. 
For  a  continuously  twice  differentiable  function  4>  on  [0,l]n,  it  is  easy  to 
verify  that  the  following  conditions  are  equivalent  (see  Lorentz  [7]): 


(3.1) 


*(u.+  h,  u.+  k)  -  $(u.+ 
i  j  i 

for  all  i*j ,  0  <  u^  ^ 


h,  u j )  -  *(ui,u>  k)  +  $(ui(  Uj)  s  0 

u.+  h  <  1,  0  s  u.  s  u.+  k  <  1. 

i  13 


(3.2) 


$(u  +  h,  u.+  h) 
i  J 

for  all  i*j , 


-  4>(u.+  h,  u.)  -  $(u.  ,u.+  h)  +  4>(u.  ,u.) 
i  3  l  3  l  3 

0  <  u.  £  u. ♦  h  s  1,  0  s  u.  s  u.+  h  <  1. 

ii  J  J 


>  0 


<3-3>  nr  3ir«“i . V20 

1  3 


for  all  i*j . 


Therefore,  due  to  the  symmetry  of  and  the  above  equivalence,  it  suffices 
to  note  that 


3u„ 


Vv 


•V 


I 

—12  ’  ^ 


U3  ...Un 

,+...♦£  =n-k 
3  n 


l-e 3  l-e_ 

(l-u3)  ...Cl-un) 


S  0 


(where  represents  an  n-2  component  vector  of  zeroes  and  ones) . 


Remark  3.2.  Let  (Xj,...,X  )  and  (Y^,...,Yn)  be  random  samples  of  size  n  from 
populations  with  life  distribution  functions  F  and  G  respectively.  Barlow  and 
Proschan  [1]  show  that  if  G  }  F  where  G  and  F  have  common  mean,  then 

(EY[1],'*''EY[n])  >  (EX[i EX[n])- 


j 


-  12  - 


Shaked  [13]  proves  the  same  result  under  the  more  general  assumption  that 
m 

G  >  F.  His  proof  uses  the  characterization  of  Corollary  2.3  together  with 
the  fact  that 


V*1 . tn)  =  1 

is  (separately)  convex  for  each  k. 


[n]+* ' *+  *[k] 

It  follows  that 


EV[n]*’"'*  EY[k]  ’  |  VS . tn)dG(t1}...dC(tnJ 

2  /  V'l . tn)dF(t1)...dF(tn) 


=  EXr  . . .+  EXr.  . 
[n]  [k] 


fa) 

Remark  3.3  Suppose  that  for  each  a  e  A,  ’  is  distribution  function 

on  R,  and  that  y  is  a  probability  measure  defined  on  a  o-field  of  subsets 

of  A.  One  may  define  the  n-variate  distribution  function (assuming  appro- 

(a) 

priate  measurability  conditions  on  F  ) 

F(x1,...,xn)  =  /  F(a)(Xl)  ...  F(a)(xn)dy(a). 

A 


If  random  variables  X^,...,Xn  have  such  a  joint  distribution  function, 
they  are  said  to  be  'positively  dependent  by  mixture*.  Given  X^,...,Xn 
positively  dependent  by  mixture,  let  Y^,...,Yn  be  independent  random 
variables  where  Y^  is  distributed  as  X^  for  i  =  l,...,n.  Shaked  [12] 

(See  also  Marshall  and  Olkin  [9]  and  Proschan  [10])  has  shown  that  in  this 
case 

(EY[1]’-**,EY[n])  >  (EX[1]...,EX[n]). 

m 

Remark  3.4  Theorem  3.1  shows  that  if  GA  >  F^^  for  all  i*l,...,n,  then  for 
n  n 

any  k  1  Yr.,  is  "more  variable”  than  T  Xr.-,  (in  the  terminology  of 
i-k  L1J  n  i=k  L1J 

Ross  [11])  or  that  £  Yr.,is  "lar8©r  in  mean  residual  life”  than  )  xr., 
i*k  L1J  i«k  ^ 


(in  the  terminology  of  Stoyan  [14]).  Since  \(t^ , . . .  ,tn)  =  t[j(] 

is  convex,  this  also  follows  by  using  the  result  that  if  Xj.-.-X  , Yj, 

. . .  ,Y  are  independent  and  Y^  is  "more  variable"  than  for  i=l,...,n, 
then  \(Yj » • • • .Yn)  is  "more  variable"  than  '*'ic(xi»*  •  •  »xn)  (see  Ross  [11]). 

Remark  3.5  If  X^,...,Xn  are  independent  HNBUE  random  variables,  then 
Theorem  3. lb  could  be  useful  in  constructing  bounds  on  the  expected  order 
statistics  EX^^, . . .  ,EXj-nj. 


Example  3.6  Let  us  consider  the  following  problem  of  general  interest. 

n  components  are  to  be  purchased  in  order  to  form  a  coherent  system  (for 

example  a  k  out  of  n  system),  and  all  of  the  components  are  to  be  purchased 

from  either  company  A  or  company  B.  Let  us  suppose  that  each  company  makes 

the  claim  that  components  of  type  i  have  mean  life  ik  (i=l,...,n),  but 

that  company  B  is  known  to  be  'more  variable'  than  company  A  in  the  production 

of  any  type  of  component.  If  we  wish  to  maximize  the  mean  life  of  the 

system,  from  which  company  should  we  buy? 

Let  X^,...,Xn  and  Yj,...,Yn  be  random  variables  representing  the 

lifetimes  of  the  components  from  A  and  B  respectively.  If  we  can  assume 

that  the  components  function  independently  within  the  system  and  that  Y^ 

m 

is  more  variable  than  X.  in  the  sense  that  G.  >  F.  (where  X.  ^  F.  and 

i  iiii 

Y^  'v  G^)  for  all  i=l,...,n,  then  we  know  that 

(EYm,...,  EY[n])  >  (EX[i;] . EX[n])- 

In  particular  EY^j-  EX[i]  *  0  81011  EY[n]  ’  EX[n]  2  °*  Therefore  our 
system  is  a  series  system  we  would  buy  from  A,  while  if  it  is  parallel  we 
would  buy  from  B.  This  result  was  observed  by  Marshall  and  Proschan  [8]. 
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For  a  more  general  k  out  of  n  system,  we  would  be  interested  in  the 
expected  order  statistics  EX^^^  and  EY^  ^+1]  order  to  compare 
companies  A  and  B.  Although 

^EYCl]***‘  >  CEX[1]»...,EX[n;i), 

EYj- ^ -j-  EXj-^-j  may  theoretically  at  least  undergo  many  sign  changes  as  i:l  -*•  n 

even  in  the  case  when  F.  =  F  and  G^=G  for  all  i=l,...,n.  However  under  the 
m 

assumption  that  G  >  F  where  G  and  F  are  continuous,  G  is  strictly  increasing 
on  its  interval  support  and  G(0)=F(0)=0,  one  may  show  that  the  number  of  sign 
changes  in  EY^-j  -  EX^  is  no  greater  than  the  number  of  sign  changes  in 
G(x)-F(x)  as  x:0  -*■»  .  Since  F*’*(t)  Fn  *  (t)  is  totally  positive  of 

order  ®  in  i  and  t,  this  follows  using  the  variation  diminishing  property 
of  totally  positive  functions  and  the  identity 

EYm  -  EXm  »  Jq  n(G_1F(t)  -  E^"1  (t)F n_i  Ct)dt 

(see  Barlow  and  Proschan  [1]).  In  particular  if  F  crosses  G  once  then  there 
exists  a  constant  C  (depending  on  n,  F  and  G)  such  that 

EYr . ,  -  EXr s  0  for  i  <  C 
Ci]  Ci] 

and 

EYj-i;J  -  EXj-.-j  2  0  for  i  >  C. 


IS 
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